Theory of vortex states in magnetic nanodisks with induced Dzyaloshinskii-Moriya 
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Vortex states in magnetic nanodisks are essentially affected by surface/interface induced 
Dzyaloshinskii-Moriya interactions. Within a micromagnetic approach we calculate the equilib- 
rium sizes and shape of the vortices as functions of magnetic field, the material and geometrical 
parameters of nanodisks. It was found that the Dzyaloshinskii-Moriya coupling can considerably 
increase sizes of vortices with "right" chirality and suppress vortices with opposite chirality. This 
allows to form a bistable system of homochiral vortices as a basic element for storage applications. 
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I. INTRODUCTION 

Physics of magnetism at nano- and submicrome- 
ter scales is an object of broad and intensive scien- 
tific investigations stimulated by various possible ap- 
plications including magnetic random access memory, 
high-density magnetic recording media, and magnetic 
sensors. [l| Nanoscale magnetic dots with vortex states 
are considered aspromising components for such spin- 
tronic devices. @, S 13 Numerous experimental observa- 
tions show that such vortices consist of a narrow core 
with a perpendicular magnetization surrounded by an ex- 
tended area with in-plane magnetization curling around 
the center (Fig. QJ. 0, H, 0, H, [§| It has been proposed 
to use both the up and down polarity, i.e. the perpen- 
dicular magnetization of a vortex, or the rotation sense 
of the curling in-plane magnetization as switchable bit 
elements in memory devices. @, H, H[ 

The vortex state in thin film elements results from the 
necessity to reduce the demagnetization energy in com- 
petition with the exchange coupling, [l^, EH In circular 
disks the axisymmetric ground state arises for diameters 
in the range of a few nanometer up to many 10 nm de- 
pending on magnetic material. For very small diameters, 
a single domain state occurs, very large film elements 
form multidomain states. Within the usual micromag- 
netic calculations the shape and size of the vortices are 
determined by the competition between exchange and 
stray field energy. [12] In particular, the vortices with dif- 
ferent chirality are degenerate: the four possible vortex 
ground-states differentiated by their handedness and po- 
larity all have the same energy within these standard 
micromagnetic models. 

The broken inversion symmetry at the surface of mag- 
netic films, however, induces chiral magnetic couplings in 
the form of Dzyaloshinskii-Moriya (DM) exchange as a 
result of relativistic spin-orbit interactions. [l3l Q, EEj] 
The chiral DM interactions destabilize collinear mag- 
netic states and are able to create a lar ge v ariety of he- 
lical and Skyrmionic spin textures. [HI Il6l IT7L I18L EH 
The mechanism and phenomenological models of the 



surface-induced DM couplings, along with possible ob- 
servable effects in magnetic films, have been discussed 
earlier. [3, EE[ These theories now are supported by 
modern quantitative ab initio calculations for magnetic 
nanostructures. [20I |2ll [H, [HI Recent experiments 0, 
|25| provide clear evidence for these surface-induced DM 
interactions, as they display long-period modulated non- 
collinear magnetic states, which can be identified as chi- 
ral Dzyaloshinskii spirals. 16] Chiral effects observed for 
magnetization processes in vortex states of magnetic nan- 
odisks may also belong to this class of phenomena. 

In this work, we describe the chiral symmetry breaking 
in the vortex ground states of circular thin film elements 
within a basic micromagnetic approach. As the vortex 
states are chiral themselves, the effect of the chiral DM 
is less obvious. However, in the presence of DM interac- 
tions the chiral degeneracy of the left- and right-handed 
vortices is lifted. The simplicity of the circular vortex 
structure makes them amenable to detailed theoretical 
investigations. Here, we calculate the differences between 
the core shapes and sizes of left- and right-handed vor- 
tices in the presence of DM couplings. These differences 
of core structure may be observable in experiments, e.g. 
as differences in core diameter or net polarity of vor- 
tices, when switching their chirality. We suggest that 
such experiments can be used to determine the magni- 
tude of surface-induced DM couplings in ultrathin mag- 
netic films/film elements. 



II. EQUATIONS AND METHODS 

The energy density of a uniaxial ferromagnet with chi- 
ral interactions can be written in the following form [l9| 
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where m is the unity vector along the magnetization 
M = M s m and M s is the saturation magnetization. 
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The couplings are given by the exchange stiffness A. The 
anisotropy axis a is taken perpendicularly to the disk sur- 
face, and K u is the anisotropy constant which must be 
positive in easy-plane materials. H is the external mag- 
netic field, and H m is the stray field. The Dzyaloshinskii- 
Moriya e nerg y wd is described by so-called Lifshitz 
invariants, [l6J energy terms linear in first spatial deriva- 
tivs of the magnetization, 
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The functional form of this energy is determined by the 
symmetry of the surface/interfaces, [lij In this paper we 
use wd terms in the following form 
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which gives the allowed Lifshitz invariants of the mag- 
netization in symmetries from Laue classes 32, 42, and 
62. This wd term favours the curling mode of the mag- 
netization, where the rotation sense is determined by 
the sign of the Dzyaloshinskii constant D. Depending 
on crystallographic symmetry, other types of Lifshitz in- 
variants may occur which may favour differently twisted 
non-collinear magnetization structures. For example, Lif- 
shitz invariants (l[ x J +L < fy), possible in Laue classes 3m, 
4m, and 6m, induce a cycloidal rotation of the magne- 
tization vector. [19j The various possibilities can be de- 
duced from the corresponding list of Lifshitz invariants 
for three-dimensional Laue classes in Ref. where also 
the structures of the corresponding Dzyaloshinskii spirals 
and circular vortex-like Skyrmions have been presented. 
For simplicity, we restrict discussion here to the form of 
Eq. ©. 

The equilibrium configurations of m are derived by 
solving the equations minimizing the energy (fTJ) together 
with equations of magnetostatics. To describe vortex 
states in a disk of radius Rd and with zero or perpen- 
dicular applied field, we consider axisymmetric distri- 
butions of the magnetization and express the magneti- 
zation vector m in terms of spherical coordinates and 
the spatial variables in cylindrical coordinates: m = 
(sin 9 cos ip; sin 9 sin ip; cos 9) , r = (p cos ip; p sin p; z) (Fig. 
[TJ b). The vortices are characterized by the magneti- 
zation direction in the center, polarity p = ±1, and 
the chirality of the magnetization structure, chirality 
c = ±1. Four different vortex states can be described 
by the indices (p; c) (Fig. [TJ c). One can introduce 
the chirality of the Dzyaloshinskii-Moriya coupling as 
D = \D\c. Then the Dzyaloshinskii-Moriya interactions 
favours states with c — c and suppresses those with op- 
posite chirality c = — c. In this paper we assume for def- 
initeness that the vortices have positive chirality, c = 1, 
and study their properties for D^O. Thus, positive DM 
couplings favour these vortices, while they are suppressed 
in the opposite case, D < 0. 

Due to the non-local character of stray-field interac- 
tions the micromagnetic problem Eq. (fTJ) constitutes a set 




FIG. 1: (a) Vortex states in a circular magnetic nanodisk 
of radius Rd with axisymmetric magnetization structure, (b) 
Geometry and definition of variables of the problem, (c) The 
four possible vortex states are characterized by the indices 
(polarity p = ±1; chirality c = ±1) Dzyaloshinskii-Moriya 
couplings with chirality c = 1 (D > 0) widen vortices with 
the same chirality, c = I , and squeeze vortices with opposite 
chirality, c = — 1. 



of integro-differential equations, llj In order to simplify 
this problem we consider the limit of a thin film where 
the magnetodipole energy has a local character and re- 
duces to a "shape" anisotropy K m = 2nM^ .[27[ This can 
be added to the uniaxial anisotropy K u yielding a redefi- 
nition of the anisotropy energy in Eq. ([TJ) by an effective 
anisotropy constant K. We also introduce the character- 
istic (exchange) length l e , the anisotropy field H a , and a 
critical value of the Dzyaloshinskii constant Dq 



Do 



l e = ^/AjK, H a = 2K/M S 



K = K u + 2ttM 1 s > 0, 



(4) 



as proper material parameters of the problem. They es- 
tablish important relations between vortex solutions and 
magnetic states in laterally infinite magnetic nanolayers. 
The anisotropy field determines the equilibrium magneti- 
zation of homogeneously magnetized layers in an applied 
perpendicular field H, 



cos9 h = H/H a . 



(5) 



The constant Do gives a threshold "strength" of the 
DM coupling in comparison with the exchange and 
anisotropy: for \D\/Do > A/tt — 1.273 the magnetization 
of a layer transforms into a modulated state. [3, [n| The 
exchange length gives a characteristic radius of the vortex 
core. Most experimentally investigated nanodisks have 
radii much larger than the exchange length, Rd S> l e - In 
this case vortices consist of a strongly localized core en- 
circled by a wide ring with a constant polar angle 9 = (9/> 
(Fig.[H(a)). 
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The variational problem for functional (JTJ) has rota- 
tionally symmetric solutions ip — <p ± 7r/2, 9 — 9{p). By 
substituting the solution for ip into Eq. |T]) and integrat- 
ing with respect to ip the vortex energy can be reduced 



to the following form E = 2tt w(p)pdp, where 
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where the magnetic field H is assumed to be perpendic- 
ular to the disk plane. 

The Euler equation for 0(p) 
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with the boundary conditions 

6(0) = 0, (d9/dp) p=Rd =g(6,R d ) 



(7) 



(8) 



yield the equilibrium vortex profiles. g(9, Rd) describes 
the anchoring effect imposed by the surface energy at the 
disk edge. For D = H = g(0, Rd) = and infinite radius, 
Eq. ([7]) is related to the differential equation introduced 
by Ginzburg and Pitaevskii in their theory of supcrfluid 
vortices in liquid helium. 28] Similar equations describe 
vortex excitations in different bosonic systems including 
Bose-Einstein condensates (see e.g. (29|) and easy-plane 
magnets. [30} 

As a first step in solving the boundary value problem 
(Eqs. ([7]). l[8"])) we consider an auxiliary Cauchy problem 
for Eq. (J7J with initial conditions 



0(0) = 0, (d9/dp) p=0 = a 



(9) 



for different values of a > 0. The trajectories 9(p; a) 
for this initial value problem with < a < 00 define 
a family of solutions parametrized by a. Any solution 
9(p) of the boundary problem (UJ, © is member of this 
family for a certain value of a. A qualitative analysis 
of possible trajectories in the phase space, d9/dp = 9 p 
vs. 9, makes it possible to single out the desired solu- 
tion among other trajectories. Typical trajectories 9 P (9) 
for the Cauchy problem are presented in the phase space 
portrait shown in Fig. O For arbitrary values of a the 
lines 9 p (9) usually end by spiraling around one of the at- 
tractors (7m, 0), n — ±f , ±2 .... But for a certain value 
a the line 9 P (9) ends in the point (0/j, 0). Thus, variation 
of parameter < a < 00 (0) allows to select the solu- 
tions of the boundary problem 0, ©. The trajectories 
0(p) have arrow-like shape. Their core sizes can be intro- 
duced in a manner commonly used for magnetic domain 
walls [TT| and Skyrmions [f9[ : as the point Rq where the 
tangent at the origin point intersects the line 9 — Oh, 



Ro = h {d9/dp)- p t Q 
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FIG. 2: (a) Typical phase space trajectories for the solutions 
of the auxiliary Cauchy problem (Eqs. ifF}. (|9jl). (b) Corre- 
sponding profiles 9(p). The shown trajectories pertains to the 
case H = 0, 9^ ~ tt/2. Due to the localized character of the 
vortex states the equilibrium solutions 8(p) of the boundary 
problem (J7J, JH} are close to the localized profiles given by 
the (dashed) separatrix lines in the phase portrait. 



The solution of the boundary value problem ([Jj , ([8]) can 
be readily found from a set of solutions of the Cauchy 
problem 9(p;a). Namely, the solution is given by the 
profile 9(p) which crosses line p = R d at the angle 
[3 = arctan( 5 (0, R d )) (Fig. [21(b)). The family of the solu- 
tions 9(p; a) of the Cauchy problem establishes the con- 
nections between the equilibrium vortex size (Rq cx 1/a), 
the disk radius R dl and the anchoring energy g(9 1 R d ). 
The trajectories 9(p; a) in Fig.[5](b) visually demonstrate 
how the vortex core size depends on the disk radius Rd 
and the boundary anchoring at the edges. This depen- 
dence should be noticeable in small disks with radii com- 
parable to the exchange length l e . In real magnetic nan- 
odisks one usually has Rd ^> l c -^ In this case the solu- 
tions of Eqs. (UJ, © are very close to separatrix lines and 
the equilibrium core retains a fixed shape and size almost 
identical to the separatrix solution a and 9h, largely in- 
dependent from the radius Rd and the anchoring energy. 
In this connection see an interesting discussion on vortex 
core sizes in Ref. (3l| . 

The independence of the core structure on boundary 
conditions allows to neglect effects imposed by the an- 
choring energy at the edges. Thus, we consider here the 
problem with free boundary conditions, g(9, R d ) = 0. 
Eqs. (J7J, ([5]) include two independent material parame- 
ters, D/Do and H/H a . For fixed values of these control 
parameters the solutions of the vortex profile have been 
derived by using the following numerical procedure. The 
Cauchy problem ([7J , §§§ was solved by the Runge-Kutta 
method. Through repeated calculations for varying val- 
ues a, the correct trajectory was searched by 'shooting' at 
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the boundary condition value ((8|) . After that the profiles 
have been improved by a relaxation calculation using a 
finite-difference method for the boundary value problem, 
for details see Ref . fl9j . 

The chirality of a non-collinear structure can be mea- 
sured from the strength of the twist or helical rotation 
of the magnetization, m • ( V x m) . For the radial vortex 
structure, the local twist is given by the expression 



d6 1 

— — I — cos sin ( 
dp p 



(11) 



The sign of this expression measures the local and heli- 
cal chirality in the structure. The comparison with Eq. 
([7]) shows that the local twist is equivalent to the local 
density of the DM energy. In particular, for < 9 < it/2 
and p > the local chirality of the helical structure is 
positive. Alternatively, the local chirality can be mea- 
sured by the p — z component of the chiral current 
= (l/(8ir 2 )e plJ \m • (<9„m x d\m). The evaluation for 
the vortices structure gives a chirality 



X =(0 p /{2irp))smt 



(12) 



Both the sign of r and x can De used to determine the 
local handedness in the vortex structure. In particular, 
a change of slope for the profile from 9 p > to 9 p < 
changes the chirality of the magnetization structure (Fig. 
[3 Inset). 



III. 



RESULTS 



In this section we present the results for the vortex core 
structures first from exact numerical calculations. Then 
an analytical ansatz is discussed which offers qualitative 
insight on the dependence and mechanism by which DM 
couplings influence the core structure of vortices. Finally 
the stability and static distortion modes of the vortex 
solutions are investigated. 
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FIG. 3: Typical solutions 8(p) for different values of D in zero 
applied field. The vortices widen when their chirality c = 1 
coincides with the chirality of the DM interaction, D > 0. In 
the opposite case, D < 0, the vortex cores shrink, and the 
local chirality x(p) changes its sign outside the vortex core 
(Inset a). The core with unfourable chirality is encircled by a 
ring with negative local polarization and favourable chirality 
(Inset b). For strong negative DM coupling (D/Dq = —1.5) 
the circulation of the projected magnetization in the plane 
changes the sense of rotation from positive in the core to neg- 
ative at the edge, 9 — > —ir/2 at Rd- 
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A. Vortex solutions and magnetization profiles 

Typical solutions of Eq. ([7]) for different values of D in 
zero field are shown in Fig. [31 The effect of applied per- 
pendicular fields on the core structure is demonstrated 
in Fig. [4] All reported results from the numerical solu- 
tion of the boundary problem ([7]) , © are for a fixed disk 
radius Rd = 30i e . As discussed, the solutions well repre- 
sent the vortex core properties for any Rd 3> l e - Vortex 
profiles 9(p) have an arrow-like shape. For D > the 
angles 9 vary monotonically from zero at the vortex axis 
to Oh ©. In this case the magnetization has everywhere 
a local chirality favoured by the Dzyaloshinskii-Moriya 
interactions. The core size widens with increasing D. 

For negative D the magnetization in the core has un- 
favourable chirality. As a result in a certain point p r 



FIG. 4: Typical solutions 8(p) for different values of the ap- 
plied field and the reduced values of the Dzyaloshinskii con- 
stant: D/Do = (solid lines), 0.5 (dashed lines), -0.5 (dash- 
dotted lines). 



with 9(p r ) > 9h the profile 9{p) goes through a maxi- 
mum, 9 > 7r/2 in zero field, the slope 9 P becomes neg- 
ative, and the magnetization structure changes its local 
chirality. After that, in the range p r < p < Rd the polar 
angle 9(p) monotonically approaches the limiting value 
9h- As the magnitude of the negative D increases these 
vortices transform into those with "reverse" rotation to 
-6 h (profile with D = -1.5 in Fig. [3j. For D < the 
vortex core consists of a narrow internal part (p < p r ) 
and the adjacent ring with a reverse magnetization rota- 
tion. 
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FIG. 5: Magnetization profiles m z (p/l c ) for different values 
of D. For D > the perpendicular magnetized central spot 
broadens with increasing D. For negative D an increasing \D\ 
widens the ring with negative m z and compresses a central 
spot with positive m z (Inset). 




H/H 



FIG. 6: The equilibrium sizes of vortex core Ro (|10[) as func- 
tions of reduced magnetic field H/H a for different values of 
D/Do- Inset shows Ro as functions of D/Do for different 
values of the applied magnetic fields. 



Peculiarities of the vortex profiles for different sign of D 
are reflected in their magnetization distribution (Fig. [5]). 
Positive D increases the width of the central spot with 
m z > (Fig. [5]). As a result the total perpendicular 
magnetization of the disk < m z >= J Q d m z (p)pdp is 
larger in vortices of the right positive chirality. Negative 
D squeezes the central magnetization core with m z > 
and widens the adjacent ring with negative perpendicu- 
lar magnetization (m z < 0) (Fig. \5[ Inset). The overall 
behavior of the vortex core sizes is summarized in Fig. [5] 
by displaying the dependence of the core radii i?o, as de- 
fined in Eq. (|10p . on the external field for different zero, 
positive, and negative values of D. This dependence is 
weak and almost linear up to large fields H/H a — > 1, 
where the transition into the homogeneously magnetized 



state takes place. The Inset of Fig. [6] shows that the de- 
pendence of Rq on D is almost linear both for zero field 
and for not too large positive and negative fields. 



B. Linear ansatz and analytical results for a vortex 



Vortex profiles with strongly localized arrow-like cores 
(Fig. [3]) can be described by a linear ansatz for the core 
and a flat part, 



= 9 h (p/R), 0<p<R, 

= e h , r < p<R d . 



(13) 



Integration of the energy functional Eq. ([7]) with the 
ansatz 6(p) (TP3"]) leads to the following expression for the 
vortex energy as a function of R, 



£{R) = ICR 2 — AlnR — cDR , 



(14) 



where K, = Trg (H)K, V = 4ir gi (H)D, A = 2irg 2 (H)A, 
and g Q (H) = (6^(1 + 2 cos 2 9 h ) - W h sin 261,, - 7 cos 2 6 h + 
8co80 h - l)/(20 2 ), g x {H) = (sin 2 6^ - 6 h e&n20 h + 
9l)/(48h), gi{H) = sin 2 Q hl and 9 h is the solution of Eq. 
(J5]). In Eq. (|14p we omit terms independent on parame- 
ter R. Minimization of energy (|14[) yields the equilibrium 
radii of the core for positive and negative D 



Ri,2 = u{H)l e 



where u(H) 



, (15) 



gi (H)/g (H), 
Particularly, 



v(H) 
at zero field 



^go(H)g 2 (H)/ gi (H). 
u(0) = 1.856, v(0) = 0.988. 

Eqs. ([14]) . (fT5|) offer an important insight into the 
physical mechanism that underlies the formation of the 
vortex states. The exchange energy of the vortex core 
does not depend on its size. This reflects ageneral prop- 
erty of vortex and 2D Skyrmionic states. [3 H3| The ex- 
change energy of the adjacent ring (oc —AhiR) favours 
the extension of the vortex cores while the magnetic 
anisotropy energy oc KR 2 tends to compress them. For 
D = the balance between these ener gy con tributions 
yields the equilibrium core sizes R oc y/ A/K, equal for 
the vortices of opposite chirality. Finite values of D vi- 
olate chiral symmetry of the solutions (fl~5} stabilizing 
vortices with different sizes of the core. The difference 
between core sizes in vortices with different chirality can 
be readily derived from Eq. (Til]]) , 



(16) 



AR=\R 1 -R 2 \=2lJ-^-u(H). 

i->0 



The numerical calculations of AR also reveal a linear 
relation AR = a(D/Do), where the coefficient a de- 
pends on the applied field. Particularly a(H/Ha) has 
the following values: a(0) = 2.4, a(0.25) = 2.66, and 
o(-0.25) = 2.24. 
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FIG. 7: The first three excitation modes for D/Do = 0.5 and 
H/H a = 0.25 and corresponding eigenvalues: Ai = 2.966, 
A2 = 8.428, A3 = 15.685. The structure is radially stable 
because the smallest eigenvalue Ai is positive. Inset shows 
the first eigenvalue as a function of the applied field and for 
different values of D/Do- 




FIG. 8: The first excitation mode £i(p) at zero field and for 
different values of the Dzyaloshinskii constant. Inset shows 
£i(p) for D — and different values of the applied field. 
Corresponding eigenvalues Xi(H/H a ): Ai(— 0.5) = 0.9556, 
Ai(0) = 1.7124, Ai(0.5) = 2.3644. 



C. Radial stability of the solutions 



with 



To study radial stability of vortices we consider a 
small arbitrary radial distortion of the solutions 

8(p) of Eqs. (O, ([8]) with the boundary conditions 
£(0) = £(R d ) = 0. We insert 6(p) = 9{p) + £(/») into 
energy functional Eq. ([7J and keep only terms up to sec- 
ond order in £(p). Because 6{p) is the solution of the 
boundary value problem the first-order term must van- 
ish, and one obtains E = E^ + E^ where E^ is the 
equilibrium energy, and 



E^ 1 



with 
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The stability problem is reduced to the solution of the 
spectral problem for functional (TP?)) (for details see Ref. 
(l9(). By expanding £(p) in a Fourier series 



£(p)=]T& fc sin(afcfl(p)), 



(19) 



k=l 



where a = ir/dh the perturbation energy E^ (P~7|) can be 
reduced to the following quadratic form 



A kj = 



[Aa 2 kj cos(ak8) cos(aj9) + 
\-G(p) s'm(ak9) sin(aj0)]pdp. 



(21) 



Radial stability of the solution is determined by the sign 
of the smallest eigenvalue Ai of matrix A: if Ai > the 
solution 9{p) is stable with respect to perturbations £(/?), 
and the solutions are radially unstable if \\ is negative. 

Numerical calculations demonstrate radial stability of 
vortex solutions for positive D in the whole range of 
the magnetic fields where these solutions exist. For 
H = 0.25, D = 0.5 the first three excitation modes 
and their eigenvalues are shown in Fig. [71 The variation 
of the first eigenfunctions and eigenvalues under the in- 
fluence of the applied field and D > is shown in Fig. 
[5] The lowest perturbation eigenfunctions are connected 
with an expansion or compression of the vortices. The 
eigenvalues and the eigenfunctions correspond to certain 
magnetic resonance modes associated with radial defor- 
mations of the vortex core. 

For D < radial stable solutions exist only below 
certain critical strength of the Dzyaloshinskii constant 
\D\ < \D cr \. For \D\ > \D cr \ vortex solutions are radi- 
ally unstable. Detailed investigations of vortex structures 
for D < beyond this threshold and their stability is be- 
yond this contribution restricted to axisymmetric simple 
vortex structures. 



IV. CONCLUSIONS 



e (2) = E A *Mi 
k,j=l 



(20) We have investigated the influence of induced 

Dzyaloshinskii- Moriya interactions on the equilibrium 
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vortex parameters in magnetic nanodisks. Both numeri- 
cal and analytical calculations demonstrate strong depen- 
dencies of the vortex structures (Fig. [3|), magnetization 
profiles (Fig. [5]), and core sizes (Fig. [5]) on the strength 
and sign of the Dzyaloshinskii-Moriya coupling. Thus, 
by switching the chirality of a vortex a change of the 
vortex profile, core size, and the perpendicular magne- 
tization takes place in the presence of a surface-induced 
Dzyaloshinskii-Moriya coupling. Existing experimental 
imaging techniques should already be sufficient to in- 
vestigate these differences in vortex states with different 
chirality. d i, 0, B i, HE IM [H, H EH The calculated 
relations between strength of the Dzyaloshinskii-Moriya 
interactions and vortex core sizes (Fig. [6] and Eqs. (p~5]> . 
(|16|) ) provide a method for experimental determination 
of the Dzyaloshinskii constant D. 

From the theoretical side our results obtained within 
a simplified micromagnetic model can be extended. The 
calculated dependence of the excitation modes in Fig. 
on the Dzyaloshinskii-Moriya interactions indicate that 
magnetic resonances of vortex cores, or more generally 



dynamical effects, may provide a route to investigate 
chiral symmetry breaking in film elements with vortex 
states. Many recent experimental [H, S, and theo- 
retical 37, 38, 39l studies have been performed on vortex- 
core dynamics and magnetization reversal in nanodisks 
under the influence of magnetic fields or current pulses. 
Effects of surface-induced Dzyaloshinskii-Moriya interac- 
tions have not been considered yet for these nanomag- 
netic processes. Our investigation on static vortex struc- 
tures is a first step towards detection of these chiral cou- 
plings in magnetic nanodisks. 
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